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ABSTRACT 

Finite  element  analysis  software  used  to  calculate  the  Stress  Intensity  Factors , KI  and  KII , for  a central  crack 
in  a plate  subjected  to  uniform  tensile  load  for  different  crack  lengths  and  orientations.  Also  for  inclined  crack  the  SIFs 
of  kinked  crack  was  investigated.  For  acceptance  of  F EM  model  used  in  computation  analysis , Numerical  results  were 
compared  with  theoretical  results  which  getting  by  solutions  of  selected  equations  and  good  agreement  had  been  found 
between  them.  The  present  study  shows  that  the  main  important  role  affects  on  stress  intensity  factors  is  the  inclination 
crack  angle  (p).  For  kinked  crack,  both  of  Mode  I & Mode  II  of  SIFs  are  strongly  depend  on  the  value  of  (p  + a)  and 
there  is  no  effect  found  when  one  of  them  (p  or  a)  change.  Furthermore  maximum  value  of  Mode  II  of  SIF  of  kinked 
crack  is  found  at  about  [(p  + a)  = (50o  - 60o)]. 
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INTRODUCTION 

The  presence  of  cracks  may  weaken  the  material  such  that  fracture  occurs  at  stress  much  less  than  the 
yield  or  ultimate  strength.  Fracture  mechanics  is  the  methodology  used  to  aid  in  selecting  materials  and  designed 
components  to  minimize  the  possibility  of  fracture  from  cracks.  Fracture  mechanics  is  based  on  the  assumption 
that  all  engineering  materials  contain  cracks  from  which  failure  starts.  Cracks  lead  to  high  stress  concentration 
near  the  crack  tip;  this  point  should  receive  particular  attention  since  it  is  here  further  crack  growth  takes  place. 

There  are  three  modes  of  loading  can  be  applied  to  a crack.  These  load  types  are  categorized  as;  Mode  I 
(opening  mode),  Mode  II  (sliding  mode,  in  - plane  shear)  and  Mode  III  (tearing  mode,  antiplane  shear).  The  most 
critical  mode  is  Mode  I because  the  crack  tip  carries  all  the  stress  whereas  in  another  two  modes  (Modes  II  and  III) 
some  of  the  stress  is  carried  by  interaction  of  the  opposing  crack  faces,  Thus  Mode  I is  the  most  common  load  type 
encountered  in  engineering  design. 

Irwin  [ 1 ] proposed  the  description  of  the  stress  field  at  a head  of  a crack  tip  by  means  of  only  one 
parameters  that  called  stress  intensity  factor , K , which  could  uniquely  define  the  stress  state  at  the  crack  tip, 
without  the  need  to  determine  the  actual  stress  components  (gxx,  oyy  oxy).  Thus; 

K = a^/na  (1) 

But  on  otherside,  Equation  (1)  is  for  the  special  case  of  idealized  crack  in  an  infinite  plate.  So  that  and 
because  of  real  cracks  are  affected  by  the  geometry  of  the  component,  the  applied  load  stress  field  and  others 
factors,  equation  (1)  can  be  generalized  as  : - 
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K = Yg  Vtt a (2) 

Where  Y is  geometrical  correction  factor  which  maintain  the  uniqueness  of  the  stress  intensity  factor  by 
accounting  for  the  particular  geometry. 

Central  cracked  finite  plate  with  and  without  kinked  is  a common  specimen  in  research  and  practice  for  fracture 
mechanics.  It  has  been  studied  by  M.A.  Meggiolaro  et.  al  [ 2]  who  used  specialized  finite  element  (FE)  and  fatigue 
assessment  software  to  explain  fatigue  crack  and  crack  growth  retardation  mechanisms.  They  obtained  the  crack  path  and 
associated  stress  intensity  factors  (SIF)  of  crack  and  kinked  cracks  for  several  crack  lengths  and  angles.  Scott  A.  Fawaz  [3] 
have  been  used  the  well  - known  Newman  / Raju,  K solution  to  predict  fatigue  life  of  unsymmetric  corner  cracks  at  a hole 
subjected  to  tension,  bending  and  bearing  load.  Also,  M.  Y.  He  et.  al  [ 4 ] studied  the  mixed  mode  intensity  factor 
distributions  at  the  edge  of  semi  - circular  and  semi  - elliptical  surface  cracks  aligned  perpendicular  to  the  surface  solid 
subjected  to  remote  shear  parallel  to  the  plane  of  crack.  Andrea  Spagnoli  et.  al  [5]  discussed  the  influence  of  the  degree  of 
crack  deflection  on  the  fatigue  behavior  by  using  a theoretical  model  of  a periodically  - kinked  crack.  By  using  the  Mellin 
transform  technique,  P.  S.  Theocaris,  G.  N.  Makrakis  [6],  studied  the  stress  intensity  factors  at  a kinked  crack  - tip.  Ameen 
Ahmed  Nassar  [7],  used  new  crack  extension  technique  to  evaluate  critical  stress  intensity  factor  (KIC)  for  cracked  plate  in 
tension  for  different  crack  configuration. 

Zheyuan  Hum  et.  al  [8]  investigated  the  effect  of  crack  orientation  and  boundary  condition  on  mixed  - mode  SIFs 
of  finite  plates  with  one  and  two  collinear  cracks  by  using  boundary  collocation  method  (BCM).  R.  K.  Bhafgat,  et.  al  [9] 
have  been  studied  the  effect  of  crack  inclination  angle  on  fracture  parameters  such  as,  stress  intensity  factor  for  modes  (I 
and  II)  problem  using  FEM.  P.  Baud,  et.  al  [10]  proposed  a scheme  to  compute  the  interaction  effects  between  two 
randomly  oriented  crack  by  using  an  iterative  technique  based  on  the  method  of  successive  approximations. 

By  using  a crack  surface  displacement  extrapolation  technique,  Nathera  Abdual  Hassan  Saleh  [11]  calculated 
Mode  II  of  stress  intensity  factor  of  several  crack  configuration  in  plate  under  uniaxial  compression.  A rectangular  plate 
with  inclined  cracks  of  different  crack  lengths  at  different  crack  inclination  angles  under  biaxial  loading  condition  are 
analyzed  by  Rahul  Kumar  Bhagat  et.  al  [12]  in  mixed  mode  condition  using  finite  element  methode  (FEM)  to  determined 
the  stress  intensity  factor. 

MATERIALS  AND  METHODS 

Figure  1 illustrates  the  geometry  of  investigated  plate  with  central  cracked.  Plate  is  made  up  of  carbon  steel  with 
young's  modulus  of  (200  GPa),  poisson's  ratio  of  (0.28)  and  density  of  (7820  Kg/  m3)  [13].  It  subjected  to  uniform  normal 
stress  (a).  Considering  homogeneous  and  isotropic  material,  square  thin  central  cracked  plate  is  used  in  finite  element 
modeling.  The  length  and  the  width  of  the  plate  are  (2h  = 2w  = 200  mm).  The  crack  length  is  (2a).  For  inclined  crack,  the 
crack  angle  is  (p)  whereas  the  kinked  angle  is  (a)  as  shown  in  Figure  2. 
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Figure  2:  Geometry  of  Inclined  Central  Cracked  Plate  with  Dimension  Plate  with  Kinked  Cracked 

To  compute  the  required  results,  PLANE  183  element  which  is  shown  in  Figure  3,  has  been  used  in  the  analysis. 
The  solid  modeling  is  used  to  generate  the  2D  model  and  mesh.  It  provides  more  accurate  results  for  mixed  (quadrilateral  - 
triangle)  automatic  meshes  and  can  tolerate  irregular  shapes  without  much  losing  of  accuracy. 


K 


Figure  3:  PLANE183  Element  Geometry,  Coordinate  System  and  Node  Locations  [14] 

Furthermore,  a deformed  shape  of  selected  Ansys’s  models  for  central  crack  with  /without  kinked  crack  at 
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different  crack  angle  (P)  and  kinked  angle  (a)  are  shown  in  Figure  4. 


Figure  4:  Deformed  Shape  of  Selected  Anysis’s  Models  for  Central  Crack 
with/without  Kinked  at  Different  Crack  Angle  (p)  and  Kinked  Angle  (a) 

THEORITICAL  CALCULATION 

Among  others,  We  selected  the  following  two  different  equations  for  theoretical  calculation  of  central  crack  and 
to  compare  their  results  with  results  have  got  by  Numerical  Method.  The  selected  Equations  are  [7,  15]:  - 

+ 0.043  (-)  - 0.491  f-)2  + 7.125  (-f  - 28.403  (-)*  + 59.583  (-f  - 65.278  (-)6  + 

Vw/  Vw/  Vw/  Vw  / Vw/  Vw/ 

Impact  Factor  (JCC):  4.7204  Index  Copernicus  Value  (ICV):  3.0 
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(3)  [15] 
(4)  [7] 


Where : Kj  is  the  Mode  I of  Stress  Intensity  factor  ( SIF ),  a is  f/ie  applied  tensile  stress , a is  half  length  of  crack, 
and  w is  /7z<?  half  width  of  plate. 

When  a crack  is  oriented  with  an  angle  (p  ) from  the  positive  x - axis  as  shown  in  Figure  2 stress  intensity  factors 
of  Mode  I & II  are  as  follows  [16]:- 


KIp  = KI(o)  cos2p 
Klip  = KI(0)  cosp  sinp 

Where : KI(o)  is  the  mode  I of  stress  intensity  factor  when  (fi  = 0°)  and  is  taken  as  constant. 


(5) 

(6) 


Whereas  for  kinked  crack  with  an  angle  (a)  from  the  original  direction  of  inclined  crack  (the  plane  of  crack)  (see 
Figure  (2),  The  local  Mode  I and  Mode  II  stress  intensity  factors  at  the  tip  of  kinked  crack  (SIFs)  are  written  as  following 
below  [16] 


KIa  = CuKI  + C12  Kn 
Kna  = C21  Ki  + C22  Kn 


(7) 

(8) 


Where:  KIa  and  KIIa  are  the  stress  intensity  factors  at  the  tip  of  the  kinked  crack  and  and  Kn  are  the  stress 
intensity  factors  for  the  main  crack  (for  present  investigation  Kj  = Kip  and  Kn  = Klip),  and 


Cn  = jCos(^)+icos(^) 


Ci2  = - l [sin(  | ) + sin(j)] 
C21  = \ [sin(  | ) + sin(j)] 


C22=icos(^)+^cos(^) 


(9) 

(10) 

(11) 

(12) 


THE  VALIDITY 


Validation  of  the  finite  element  approach  with  results  available  in  the  literatures,  experimental  results  or 
theoretical  calculation  results  which  can  get  by  solutions  of  equations  are  the  most  important  for  acceptance  of  the  FEM 
model  used  in  the  computation  [12].  In  the  present  investigations,  FEM  results  are  compared  with  results  getting  by 
solution  of  equations  (3  & 4).  From  Figures  5 and  6,  it  is  clear  that  the  theoretical  calculations  of  stress  intensity  factors,  Kj 
& Kn  of  central  crack  with  or  / without  inclination  are  very  close  to  results  of  analytical  solutions  at  different  parameters 
such  as  (crack  length  (a  / w),  applied  stress  (a),  crack  angle  (p).  A good  agreement  is  found  between  the  two  results 
especially  with  equation  (3).  This  means,  FE  modeling  using  ANSYS  software  can  be  used  to  determine  the  stress  intensity 
factors  (Kb  Kn)  for  any  crack  configurations. 


www.tivrc.ors 


editor  @tjprc.  org 


12 


Lattif  Shekher  Jabur 


2500 


^ 2000 
^1500 

■<  1000 

I' 

500 


Theo.  Equation  No.  4 
*Theo.  Equation  No.  3 
Numerical . Any  sis 


100  MPa 


0 \ 1 1 T T T T T T T 1 

0 0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 

Crack  Length  ( a -tv  ) 


A pplied  stress  (MPa) 


0 0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 


Crack  1 eugib  ( a tv  ) 

( a tv  = 0.3 ) 


Applied  Stress  (MPa  ) 


Figure  5:  Variation  of  Theoretical  and  Numerical  SIF,  Ki,  of  Central  Crack 
versus  Crack  Length  (a  / w)  and  Applied  Stress 


Figure  6.  Variation  of  Theoretical  and  Numerical  SIFs,  Kx  & Kn,  of  Inclined  Central  Crack 
versus  Crack  Angle  (P)  at  Different  Crack  Length  (a  / w). 
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RESULTS  AND  DISCUSSIONS 

Central  Crack  Without  Kinked 

Effect  of  Applied  Stress,  Crack  Length  (a  / W)  and  Crack  Angle  (P)  on  SIFs  (Kr  & Kn) 


As  expected,  Figure  7 shows,  for  non  - inclined  central  crack  (p  = 0°),  Mode  I of  stress  intensity  factor  (Kj) 
increase  with  an  increasing  of  applied  stress  and  crack  length  (a  / w).  When  the  crack  inclined  with  an  angle  (p),  for  all 
crack  length  have  been  investigated,  stress  intensity  factor,  Kb  decreases  with  the  increasing  of  crack  angle  (P)  until 
reaches  to  minimum  value  (Kj  = 0)  at  (P  = 90°),  while  Mode  II,  Kn  increases  and  reaches  it’s  peak  values  at  about  (p  = 45° 
) then  decreases  with  an  increasing  of  crack  angle  (p)  as  shown.  Also,  it’s  clearly  seen  that,  the  lower  the  crack  length  (a/w) 
the  lower  of  SIFs  (Kj  & Kn).  This  observations  are  agree  with  the  results  of  [9,  12]. 
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Figure  7:  Variation  of  Numerical  SIFs,  Ki  & Kn,  versus  Crack  Angle  (P)  at  Different  (a  / w) 

The  Relationship  between  (Kr  & Kn)  of  Central  Crack 

By  comparing  the  values  of  Kj  & Kn  at  different  crack  angle  (P)  as  shown  in  Figures  8,  at  p < 45°  the  values  of 
stress  intensity  factors,  Kb  are  bigger  than  Kn’s  values  for  all  crack  length  (a/w)  and  this  differences  in  value  decrease  with 
an  increasing  of  crack  angle  (p)  even  reaches  zero  at  (P  = 45°),  whereas  When  p > 45°  the  relation  between  Kj  & Kn  values 
is  reversed  (Kn's  values  is  bigger  than  Kj's  values).  It  means,  at  P < 45°,  the  effect  of  Mode  I,  Kj  is  the  dominant,  while 
when  p > 45°  effect  of  Kn  is  dominated.  Therefore  both  of  Kj  & Kn  have  to  take  into  account  in  design  according  to  the 
expected  inclination  of  the  crack.  This  result  is  supplement  with  results  of  [17]. 
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Figure  8:  The  Relationship  between  Mode  I & Mode  II  of  Stress  Intensity  Factors  at 
Different  Crack  Angle  (P)  and  Crack  Length  (a  / w) 

Central  Crack  with  Kinked 

Effect  of  crack  angle  (fi)  and  kinked  crack  angle  (a)  on  SIFs  (KI  (tt)  & KIt  (tt)) 

From  Figure  9,  Mode  I of  SIF  of  kinked  crack,  Kj  ( a ),  decreases  with  the  increasing  of  both  crack  angle  (p)  and 
kinked  crack  angle  (a)  for  all  crack  length  (a/w),  and  the  value  of  Kx  ( a } equal  to  zero  when  (P  + a)  about  equal  to  (85°  - 
90°).  It  is  clear  that,  for  all  crack  length,  when  crack  angle  (P)  increases,  SIF,  Kj  (a}  reaches  zero  value  at  lower  kinked 
crack  angle  ( a ),  So  that,  the  effect  of  crack  angle  (p)  and  kinked  crack  angle  (a)  on  SIF,  Kj  ( a } can  be  summarized  as  : 
For  any  value  of  («),  If  crack  angle  (/?)  increases , kinked  crack  angle  (a)  decreases . It  is  an  inverse  relationship 
between  crack  (p)  and  kinked  crack  angle  (a). 
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Figure  9:  Variation  of  Mode  I of  Stress  Intensity  Factors,  Kj  (a ),  versus  Kinked  Crack  Angle  (a) 
at  Different  crack  Angle  (P)  and  Crack  Length  (a  / w). 

Kn  ( a ) takes  the  same  behaviour  of  Mode  II,  SIF  of  the  main  crack,  Kn  ( p ),  where  increases  with  increasing  both 
of  crack  and  kinked  angles  (p  & a)  until  reaches  it's  maximum  value  at  (P  + a)  equal  about  (55°  - 60°)  then  decreases  as 
shown  in  Figure  10.  It  means,  the  maximum  value  of  Mode  II  of  SIF  of  kinked  crack  Kn  ( a ) doesn’t  occur  at  (a  = 45°)  as 
had  been  found  for  main  crack,  KII(  (see  figures  6 & 7),  but  also  it  depends  on  the  value  of  crack  angle  (p).  It  is  strongly 
depends  on  the  summation  of  (P  & a). 
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In  addition  to  that,  the  main  interesting  observation  which  can  be  seen  is  ; at  (p  + a)  < (55°  - 60°),  the  higher 
value  of  crack  angle  (p),  the  higher  value  of  Kn  ( a } is,  while  at  (p  + a)  > (55°  - 60°),  the  Kn  ( a } value  reverses  where  the 
higher  value  of  crack  angle  (p),  the  lower  value  of  Kn  (0i)  is. 
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Figure  10:  Variation  of  Mode  II  of  Stress  Intensity  Factors,  KII(a  ),  versus  Kinked  Crack  Angle  (a) 
at  different  Crack  Angle  (P)  and  Crack  length  (a  / w) 

On  other  hand,  at  constant  value  of  (P  + a),  there  is  no  effect  have  been  found  in  values  of  Kj  ( a ) and  Kn  ( a ) by 
varying  the  values  of  one  of  (P  or  a)  angles,  while  there  is  a considerable  effect  can  be  seen  when  (p  + a)  change  together 
as  shown  in  figure  1 1 . 
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Figure  11:  Effect  of  Change  Crack  Angle  ( P ) or  Kinked  Angle  (a) 
on  Stress  Intensity  Factors,  Kx  ( a } & Kn  ( « ) at  Constant  (P  + a) 
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The  Relationship  between  SIFs  (Kj  ( a ) & Kn  ( a ) ) of  Kinked  Crack 

Figure  12  clearly  shows  the  relationship  between  Mode  I & Mode  II  of  SIFs  of  kinked  crack,  Kj  ( a } & Kn  ( a ). 
The  same  relationship  between  Kj  (p}  & Kn  (p>  is  found.  It  can  be  seen,  ;at(p  + a)<(55°-  60° ),  the  SIF's  value  of  Kj  ( a ) 
is  bigger  than  the  SIF's  value  of  Kn  ( a whereas  when  ( p + a ) exceeds  ( 55°  - 60° ),  the  SIF's  value  of  Kn  ( a ) is  bigger. 
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Figure  12:  The  Relationship  between  Mode  I & Mode  II  of  Stress  Intensity  Factors 
of  Kinked  Rack  at  Crack  Angle  (P)  and  Crack  Length  (a  / w) 


Finally,  figure  13  illustrated  clearly  Von  - mises  stresses  of  some  of  selected  investigated  Ansys’s  models  for 
central  crack  with  and  without  inclination  at  different  crack  and  kinked  angles  (P  & a).  From  figures  4 & 13  It  can  see  the 
effect  of  crack  angle  (P)  and  kinked  angle  (a)  on  the  stress  distribution  and  structure  deformation  in  investigated  plate 
model. 


2a  = 80  mm,  p = 0° 


2a  = 80  mm,  p = 30° 
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2a  80  mm  = , a = 30°,  p = 30° 


2a  = 80  mm,  p = 45°,  a = 45° 


Figure  13:  Von-Mises  Stresses  of  Selected  ANYSIS’S  Models  for  Central  Crack 
with  / without  Kinked  at  different  Crack  Angle  (P)  and  Kinked  Angle  (a). 

CONCLUSIONS 

The  results  of  present  investigation  can  be  summarized  as 

• The  main  important  role  that  affects  on  Mode  I and  Mode  II  of  stress  intensity  factor  is  the  crack  angle  (p)  Where  at 

P < 45°,  the  value  of  Mode  I (Kj  ( p >)  is  bigger  than  the  value  of  Mode  II  (Kn  ( p >),  whereas  when  p > 45o,  the 
value  of  Mode  II  of  SIF  ( Kn  ( p ) is  bigger  than  the  value  of  Mode  I of  SIF  ( Kx  ( p).  Furthermore  the  maximum 
value  of  Mode  I (Kj  ( p >)  is  at  p = 0°  while  the  maximum  value  of  Mode  II  (Kj  ( p })  at  (p  = 45°) 

• Both  of  Modes  I & II  of  SIF  of  kinked  crack,  (Kj  (a)  & Kj  (a)),  are  strongly  depend  on  the  value  of  (P  + a),  and 
there  is  no  effect  found  when  one  of  them  (p  or  a)  change. 

• Maximum  value  of  Modes  II  of  SIF  of  kinked  crack  is  found  at  about  [(p  + a)  = (55°  - 60°)],  but  at  (p  + a)  < 
(55o-60o),  the  lower  value  of  crack  angle  (P),  the  lower  value  of  Kn  (a)  is,  whereas  at  (P  + a)  > (55o  - 60o),  the 
lower  value  of  crack  angle  (p),  the  higher  value  of  (Kn  (a>)  is. 
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